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I» INTRODUCTION 


The dynamical state of the Ionosphere la distinctive in three 
respects! tidal oscillations achieve maximum amplitude and exceed 
the prevailing winds; Irregular winds, subject to rapid variation with 
height, attain an equal intensity; and the turbulent regime of lower 
levels gives way to a non— turbulent domain above* Extensive observational 
data of winds and tides at the ionospheric heights using meteor and rocket 
vapor trail techniques are available in the literature. Analyses of these 
data reveal strong seasonal variations In the tidal winds, both diurnal 
and semidiurnal at height near 90 km. Atmospheric oscillations, with 
periods usually no greater than an hour, are generally called atmospheric, 
gravity waves. These waves (or oscillations) are presumably responsible 
lor perturbations in the neutral and Ionized structures of the atmosphere 
both in the mesosphere and in the lower thermosphere of the terrestrial 
atmosphere. 

Hines (1960) has pointed out that the Irregular fluctuations, whose 
periods are in the 10-30 minute range seen on wind and temperature profiles 

In the 80-300 km region, are due to atmospheric gravity waves. Perturbations 
of varying intensities and spatial dimensions are observed in the electron 
density distribution of the Ionosphere. Thitheridge (1963), from the total 
electron content measurements^ showed that the moving perturbations in the 
ionosphere range from 50-500 km in size. These moving perturbations in the 





electron density distribution of the ionosphere are broadly classified as 
travelling ionospheric disturbances (TID’s). The TlD*s frequently appear 
in the form of wave-like motions or complexities on the ionosonde records. 
Several source mechanisms * namely « auroral disturbances, auroral electro- 
jets, severe magnetic storms, and jet streams have been suggested for the 
generation of TID’s. Several local weather disturbances have been suggested 
as sources for high frequency fluctuations in the F region levels. Recent 
investigations by Georges (1967), and Baker and Davies (1969), have shown 
an apparent correlation between tropospheric weather effects and coupling 
at ionospheric heights. Georges (1967) suggested that observations of a 
relatively narrow spectrum of waves with periods near 3 minutes were results 
of the atmospheric filtering of acoustic-gravity waves. Further investigations 
by Georges (1968), and Davies and Jones (1972), have shown that the follow- 
ing criteria should be satisfied in order to have coupling between the 
troposphere and the ionosphere ; 

(1) Cloud tops must be greater than 12 km. 

(2) Thunderstorm cells must lie within about 200-250 km from 
the point on the ground directly below the ionospheric 
reflection point. 

Earlier investigations of ionospheric disturbances have shown an 
apparent correlation between the neutral atmospheric perturbations and 
wave-like phenomenon observed on several ground baaed ionospheric sounding 
records. The CW (continuous wave) Doppler sounding technique developed by 
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Davies (1962) , and Davies and Baker (1966) was one of the very successful 
techniques used to monitor ionospheric disturbances on a 
routine basis. This technique was used in our present investigation of 
the atmospheric dynamics. 

In this report t we present the details of a three->dlmenslonal CW 
Doppler sounding system currently under operation at the NASA-Marshall 
Space Flight Center, Alabama, followed by a discussion on the properties 
of the neutral atmosphere and on the theory of doppler sounding technique. 
Finally, methods of data analyses used to investigate the dynamical 
phenomena at the Ionospheric heights and suggestions for future investi- 
gations are provided at the end of the report. 
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II. PROPERTIES OF THE NEUTRAL ATMOSPHERE 


The primary properties of the atmosphere are Its density, pressure, 
temperature, composition and motion. The density of the terrestrial 
atmosphere decreases exponentially with Increase in altitude. The gravi- 
tational field Is mainly responsible for making the atmosphere Inhomogeneous and 
also for Introducing anisotropy Into the propagation of atmospheric waves. 

A vertically displaced fluid parcel of air experiences a buoyant force 
which tends to restore the original equilibrium. The rapidity with which 
the equilibrium Is restored characterizes the atmospheric oscillation 
frequency. A comprehensive theory of atmospheric gravity waves at the Ionospheric 
heights In an Isothermal atmosphere was given by Hines (1960, 1968). Recently 
Wu, et al. (1974) examined the effects of wind shear on acoustic-gravity wave 
propagation In the thermosphere. When discussing the propagation characteristics 
of the acoustic-gravity waves in the thermosphere, one can neglect the effects 
due to Che earth’s rotation (the Coriolis force), since the wave periods under 
discussion are less chan two hours, and the effects of the electromagnetic field 
on large scale motions of the atmosphere. The governing equations may be 
written as (Wu, et al. , 1974) 

Continuity, 

if + V • (p y) - 0 

3t 


Momentum, 
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Energy, 
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Equation of State, 

P " PRT, 

where, p, v, T, p, g, q, ir, R, and Cy denote the density, velocity, temperature, 
pressure, gravitational acceleration, heat flux, stress tensor, universal gas 
constant and specific heat at constant volume, respectively. The heat flux and the 
stress tensor can be written as 
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with K, M and ^ being the therrael conductivity, and the first and second viscosity 
coefficients respectively. 

Choosing the Z-axls along the vertical height, the x-axis along the eastward 
direction, and the y*-axls along the northward direction, Wu et al. , (1974) 
obtained the following dispersion relation for acoustic-gravity waves. 


S'* - [c* + kg*) + 1 Y g kg] + 


(Ic^k k - <Y-1) g k ] w + (Y-l)k ^ g^ - 0 

* * dz * 
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Introducing 


K, - •'l - 1 i^f- • 

we have the following dispersion relation including the effects of wind 
velocity and wind shear 

S'* - (k,* + kz^) S* + [ic* k^ k^- (1 Y-*l)g ^ 

2 dz 

2 2 

+ (Y-1) kx^ g* - f” - 0 

4c 

When the wind velocity and wind shear vanish the above equation becomes identical 
to the result obtained by Hines (1960) . 

Hooke (1968) investigated the effects of the passage of internal 
atmospheric gravity waves on F-reglon ionization and concluded that the 
gravity wave induced changes in the values of recombination and loss 
coefficients in the F-region appear to be small. Following Hooke (1968) 

In the interpretation of CW doppler array data in terms of internal atrno^ 
spheric gravity waves, we have assumed that the mechanism for the inter- 
action between atmospheric waves and ionization preserves such wave 
characteristics as period, phase and group velocities of the wave front. 

Very large TID*a which are, in most cases, accompanied by magnetic storms 
are Interpreted by Georges (1968) as free internal gravity waves with 
properties determined by the thermospheric parameters. Estimated theoretical 
wave front tilts (Toltoy, 1963) agree roughly with the observational results 
of Georges (1968) and Thome (1964, 1968). Observational results on the 
medium scale TID characteristics appear to fit Friedman's (1966) explanation 
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In terms of Imperfectly ducted internal atmospheric gravity waves In a 
leaky duct situation and showed that, for a given mode, phase and group 
velocities are nearly equal and relatively independent of period. The 
observed medium scale TID speeds of about 300 m sec"^ and the apparent 
direction of TID's travelling away from the winter poles agree with 
Friedman's (1966) calculations. The observed seasonal variation In the 
TID propagation direction may be explained by the seasonal variations In 
the mesospheric temperatures. It Is observed that the winter polar 
mesospheric heatings are frequent and that they affect the mesospheric 
temperature height profile (Maeda and Young, 1966), The temperature profiles 
of Champion (1957) do Indicate an appreciable weakening of the winter duct 
around 80 km, compared to the summer, for the middle latitudes. These 
temperature variations In the mososphere cause leaky duct conditions which 
allow atmospheric gravity wave energy to penetrate to F-reglon levels to 
support the wave propagation. 
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m. RADIO SOUNDING OF THE IONOSPHERE 


A model C-4 type conventional swept frequency lonosonde, a pulsed 
phase path sounder, a three dimensional CW doppler sounder array, and an 
alrglow observatory are currently In operation. A summary description of 
the phase path sounding system and measurements follows: 

The electric field intensity E(C) of a radlowave at some point 5 
on the ray can be expressed as 


E(0 


A(0 Exp ( j (ojt - 


J 


y ds) ] 


where the integral is taken along the ray from some arbitrarily chosen 
origin ^ - 0. Then, at any instant t, the difference in phase A(|) between 
the signal at the point K and the signal at the origin is given by 


- k. 
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Since A<fi/2Tr is the number of vacuum wavelengths represented by the difference 
in phase, we see that the integral 


P 


J 


y ds 


o 

yields a distance equal to the total length of the specified number of vacuum 
wavelengths. This quantity p is called the phase path of a radio wave 
passing through the ionosphere. 

Two different techniques are presently employed for making phase 
path measurements. The first method consists of transmitting fixed 
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frequency pulses with coherent detection used to indicate the phase of the 
echo relative to a stable reference. This technique also provides accurate 
determination of the group paths. Measurements of phase paths and group 
paths of radio waves as a function of frequency are presently being made 
using this phase path sounder system. The pertinent specifications of the 
phase path sounder system are listed in Table 1. 

The second method of measuring phase is by using a continuous wave 
doppler sounding technique. When a radio wave of a constant frequency f 
is transmitted from a stable transmitter and received via ionospheric 
reflection, it is observed that, the frequency of the received signal 
fluctuates slightly. The frequency fluctuations for ionospheric propagation 
are usually of the order of a few Hertzs. These perturbations are related to 
changes In the phase path of the probing signal. The relation connecting 
frequency fluctuations and phase path variations in the absence of earth's 
magnetic field can be expressed, as (Davies, 1965) 


where 


Af 


-f dp 
c dt 


f - transmitted frequency, c « speed of light in free space, 

P “ pds “ phase path 

and U Is the refractive index along the ray path s. The phase path, and hence 
A£,are functions of both p and a. Several relations connecting the doppler 
frequency variations and ionospheric parameters are summarized in Table 2. 
Contributions to the frequency variations may occur along the radio path 
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within the Ionosphere, but due to the 1/y factor, the region near reflection 
seems to be the most effective for producing the perturbations In the probing 
frequency. A frequency shift can result from either of the two following 
effects, or a combination of both. 

(1) change In ray path (s) caused by change In the height of 
reflection, or 

(2) change in refractive Index (y) along the path caused by 
change In electron density below the level of reflection. 

Considering the first case. If we assume mlrror-like reflection at a height 
h, with no changes In refractive Index below h, then 

. -2 — 

dt dt 

Af 2f dh 

Af « — — 
c dt 

for a given rate of change of reflection height 

A f a f 

In the second case. If we assume a constant height of reflection and change 
In electron density below the level of reflection then, 

Afa 1/f 

In reality the situation Is more complicated, but the Important point Is 
that the variation of Af with f gives an Indication of the cause of the 
doppler shift Af. An Important case where doppler shifts are not dominated 
by reflection heights Is that of solar flare Ionization enhancements due 
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Che radiowave and the atmospheric disturbance Is not known. This can be 
resolved from lonosonde measurements. 

Davies, et al. , (1969) showed the analogies between some magneto- 
ionic ray formulas and corresponding formulas for neutral atmospheric waves 
in an Isothermal atmosphere. The acoustic waves are analogous to radio 
waves of frequencies greater than the electron gyrofrequency whereas 

atmospheric gravity waves are somewhat similar to whistler waves. Some 
of the comparisons between atmospheric waves and radio waves are presented 
in Tables 3 and 4. 
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IV. DATA OVERVIEW AND METHODS OF ANALYSIS 

The visible record of the doppler data la obtained by frequency 
analysis of the signal. The frequency analysis Is carried out by using 
a UA6B Ubiquitous audio spectrum analyzer. The output from the spectrum 
analyzer Is displaced on a facsimile chart. About two weeks of data tape 
may be played back producing about 20 feet of facsimile chart record. It 
Is not proposed to analyze the complete data that are being collected, but 
Instead to study specific atmospheric- Ionospheric events seen clearly on 
the doppler record. The initial 20 feet of facsimile chart record con- 
taining approximately two weeks of data are being used to Identify specific 
events such as travelling Ionospheric disturbances, and thunderstorm, solar 
flare and magnetic storm effects. After visual identification of the 
events, the segments of the data tape containing the events are played back 
at a dlower speed for better time resolution. By playing back the tapes at 
a slower speed and recording the data, one can achieve a time resolution of 
about five seconds to study the fine structure In the doppler variation. 

The specific Identified segments of the analog doppler data are recorded 
on a separate tape and sent to the computer laboratory of NASA-MSFC. 

These data are then digitized at ten second intervals for further analysis. 

The geographic location of the array Is shown in Table 5 . The spaced 
transmitters are located at Muscle Shoals, Alabama, at Nickajack Dam, Tennessee, 
and at Fort McClellan, Alabama. As mentioned earlier, one of the alms of the 
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program is to determine the horizontal and vertical trace velocities of 
the TIDs and then determine the speed, elevation and azimuth of a 
travelling ionospheric disturbance from the recordings of all nine 
field transmitters. The following methods are used to derive the 
characteristics of the travelling ionospheric disturbances. 

The horizontal and vertical phase velocities are determined from the 
doppler records by determining the relative time displacements. Imagine that a 
plane wave front is moving across the array system as shown in Figure 1. 

Let Tjjjj and be the time delays between the reflection points M-N and 
M-F respectively. Here M, N, and F refer to the ionospheric reflection 
points corresponding to the transmitters located at Muscle Shoals, Alabama, 
Nlckajack Dam, Tennessesk and Fort McClellan, Alabama, respectively. The 
apparent trace velocities and Vj^p are given by the following relations: 

MN MF 

" Tmn “ TMF 

The magnitude of Vjj and the direction (j> (measured east from north) can be 
determined from the following set of equations 


~ ^MF ^MF 
Vmn Cos - Vj^p Cos ii>^p 


and the magnitude of Vji is given by 
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where end ere the engles between the geogrephlc north and the lines 
MN and MF respectively. The true velocity siay be determined by computing 
the vertical trace velocity Vy. The vertical trace velocity Vy may be 
determined from the time displacements at various frequencies and the 
difference of reflection points for a given transmitter and receiver 
location* are assumed to lie on a vertical plane* The heights of reflections 
for various frequencies Involved are determined from the lonogram analysis. 
From the values of horizontal phase velocity and vertical trace velocity 
Vy* the true phase velocity may be determined as shown In Figure 2. 

For a given sample of the doppler sounder record, the time displacements 
are determined by computing the cross-correlation functions. Consider two 
time series functions given by f(t) and g(t). The cross-correlation 
coefficients of f(t) and g(t) are given by 


Cfg (T) 

f and g are the time averages of f(t) and g(t). 

T Is the time displacement between the two functions being correlated. 
Sf and Sg are the standard deviations of f (t) and g(t) where 


^f ^g 




J 


lf(t) - f] [g(t 


+ T)-g] dt| 


1/2 


^ f 

r^o J 


[f(t) - fi^ dt 


to 


14 



1/2 


’8 


‘ 

‘'l ’ ’'o J 


Ig(t) - g] dt 


The auto-'correlatlon coefficient Is given by 


M - ^ f [f(t> - f) [f(t+T)-fI dti 

(s,)’ I ti- to J I 

The time interval t^^ la the period the Ionospheric disturbance lasted 

on the doppler sounder record. 

The characteristics of some of the wavelike ionospheric disturbances 
seen on CW doppler sounder records are subjected to power spectral analysis 
to Investigate the wave characteristics. One of the characteristics of the 
Ionosphere Is the electron density variability with time at a given place. 

The Ionosphere may be treated as a quasi-statlonary random phenomenon since 
several unknown processes are In operation In any given situation. Power 
spectral analysis Is often used to Investigate the characteristics of the 
quasl~statlonary random process. Tolstoy and Montes (1971) successfully used 
power spectral analysis of CW doppler sounder data to Investigate the phase 
height fluctuations in the Ionosphere between 130 and 150 km. The general 
problon of spectrum analysis was formulated by Blackman and Tukey (1959) . 

Consider a cross-correlation function of two time series 


4T/2 

fl(t) • dt 

-T/2 
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where U) Is the angular frequency In radians per unit time. In general the 
cross correlation function Is asynunetrlcal and therefore the cross spectrum 
function would be complex 

P 12 ^ ^ ^uad 

where Cq(<j)) Is the real part and Quad (o>) Is the imaginary part of the 
cross spectrum function. 

The phase relationship between frequencies Is given by 
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IQ j (ui) 

A 0 (w) « 0_(w) - 0, (oj) ■ tan“^ 

2 1 L Co (w) J 

where 6 j^(to) phase of time series 

62 (<i)) phase of time series f 2 (c) 

The coherence estimates are given by 

Pjl(«) • P22<‘^) 

where |y( 1 u)| is the magnitude of the complex coherence. 

The amplitude of the cross spectrum at each frequency Is the product of 
the corresponding amplitudes in the two time series. If a particular spectral 
component is absent from either series, then It will be absent In the cross- 
spectrum and the coherence will be a minimum at that frequency. 

Blackman and Tukey (1959) have recommended that In order to minimize 
leakage from the low frequencies to the high frequencies, the data should 
be pre-whltened prior to the power spectrum analysis. This may be accomplished 
by applying to the data a simple derivative filter. The smoothed spectral 
estimates may be obtained by using a 'banning window' which has side lobes 
of the order of 1% of the main lobe. The number of lags are usually equal 
to 10 % of the data sample giving 20 degrees of freedom. 

The lowest frequency that can be analyzed in the data sample is 
determined by the maximum lag, since f^tin “ (2T)“^. The digitization 

Interval, At, determines the highest frequency, fmax “ (2At)“*. The pre- 
whltenlng filter may be restored by multiplying each frequency estimate 
by the Inverse of the square of the amplitude response at that frequency. 
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V. ATMOSPHERIC-IONOSPHERIC DISTURBANCES 


A summary of the observations using NASA-MSPC Doppler Sounder array 
Is presented and discussed here. The experimental results obtained from 
New York Doppler Sounder array are also analyzed and the results are 
presented here. The technical details and geographic location of the New 
York Doppler Sounder array are given In a report by Rao (1972). 

A typical doppler sounder recorder from the NASA/KSFC array , showing 
the doppler variations (phase path variations) over a 24 hour period for 
the operating frequency 4.0125 MHz Is shown In Figure 3 . It can be seen 
from the figure that during the day time the doppler variations are quite 
narrow with small and rapid frequency fluctuations. Sunrise effects on the 
doppler record are shown around 0530-0600 CST. During sunrise conditions, 
there Is an Increase In the electron density of the Ionosphere due to 
photo -Ionization. The Increase In the electron density causes a decrease 
In the radio phase path length, which appears as an increase In doppler 
shift in frequency. After the sunrise effects, the doppler variations remained 
fairly constant with localized small rapid frequency fluctuations superposed 
on the record. These rapid frequency fluctuations are mainly due to the 
changes in the electron density near the reflection levels caused by the 
rapid production and loss mechanism of electrons due to the solar radiation 
and the movement of localized small scale Ionospheric irregularities. The 
doppler sounding technique of the Ionosphere Is capable of separating the 
magneto -lonXc componentsof the electromagnetic wave propagating in the 
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Ionosphere when there Is a large doppler shift in the observations. The 
separation of themagneto-’lonic components is seen from Figure 3 in which 
a weak trace is shown on the main doppler trace during the periods between 
1800-0400 CSU indicating the presence of an extraordinary magneto-ionic 
mode. It can also be seen from the figure that the doppler variation 
between 1800 and 0200 CST shows large fluctuations indicating the presence of 
wave— like phenomena. Such events seen on the doppler records are classified 
as travelling ionospheric disturbances (TID's). Examination of several 
days of doppler records taken at fixed frequency transmission under quiet 
ionospheric conditions indicate a general trend that the variations are 
small during day time hours whereas the night time variations are large. 

This variability in the doppler fluctuations may be due to the reflection 
height changes from the daytime to the night time on fixed frequency 
sounding records and the predominant occurrence of TIDs during sunset and 
night time conditions. On the high-frequency soundings (especially when 
operated near the F2 layer critical frequency) the doppler signal is either 
lost or a diffuse trace is recorded. Under spread F conditions, the doppler 
record becomes very diffuse and a wave like pattern is noticed in the 
diffused (often lost) trace of the doppler record. These investigations are 
currently under progress from the data that are being collected using the 
NASA/MSFC doppler sounder array. 

Travelling Ionospheric Disturbances 

Some of the typical TIDs observed on our doppler sounder array are 
shown in Figures 4, 5, 6, 7 and 8. Usually the TIl/s are observed as a wave like 
fluctuation on the doppler sounder record. The disturbance first appears on 
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the high frequency sounding records followed by middle low frequency 
sounding records. It Is also observed from the typical TID's recorded 
on the doppler system that the amplitude of the doppler fluctuations 
are larger for higher sounding frequencies and lower for the lower 
sounding frequencies. The TIDs are generally interpreted as manifestations 
of Internal atmospheric gravity waves (Hines, 1960). The characteristics 
of the Internal gravity waves are that the amplitude of the waves Increases 
with altitude. This Is observed on the doppler records as the perturbation 
amplitude Increased with the probing frequency. The CW doppler data for 
several TID events for the MASA/MSFC array are under analysis and the results 
will be reported later. 

During the period of this Investigation, we have analyzed some of 
the TID data collected from a CM doppler sounder array located In the New 
York area (41°N, 74°W) and the results are reported here. The horizontal 
phase velocities for the TID events observed during winter (November, 
December 1969 and January 1970) and equinoctial (October 1969, February, 
April 1970) conditions have been computed. The diurnal variation of phase 
speeds for winter and equinoctial seasons are shown in Figures 9 and 10 > 
respectively. The daytime winter speeds range from 70 to 140 m sec~^ 
whereas the Aight time winter values are higher and In the range of 140- 
340 m see”’. There are not many day time equinoctial TID events available; 
but » In generalithe day time speeds are lower than the night time. Also 
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comparison of figures Indicates that TID activity Is more predominant In the 
winter than at the equinoxes. TID's were observed on the CW doppler records 
from sunrise to late midnight during the winter whereas during the equinoctial 
months, the major TID activity occurred during the sunset and post sunset 
periods. The diurnal variation of TID propagation (directions heading) 
during the winter and equinox seasons are shown in Figures 11 and 12 
respectively. In winter TID's propagate to the Southeast during daytime and 
towards the Southwest during the nighttime. It appears, In general, that the 
TID's propagate generally in the north-south direction during winter seasons. 
The predominant propagation direction Is between west and southwest during 
equinoctial conditions. 

TID observations In the northern and southern hemispheres, as reported 
by various Investigators, are sumnarlzed In Table 6. Munro (1950, 1953), 
from an extensive study of medium scale TID’s for southern hemisphere 
locations, reported that the mean directions of travel were northward In winter 
and eastward In the equinoctial months. The horizontal speeds varied from 
85 ~ 170 m sec"^. From the results of the present Investigation, as well as 
those of earlier Investigations, It appears that during winter and summer 
periods TID's travel away from the winter poles. During equinoctial seasons, 
the TID's In the northern hemisphere have an apparent westwood motion while 
those In the southern hemisphere have an eastward motion. Examinations of 
Table 6 also show that large scale disturbances are either accompanied by 
a magnetic storm or associated with some form of auroral zone activity. On 
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the other hand* medium scale TID's are not generally correlated with the 
auroral or magnetic activity. 

Severe magnetic atoms and auroral electro jets were shown to be 
sources for large scale TIDs, (Hunsucker and Treten, 1967; Flock and Hunsucker, 
1968; Thome, 1968). The results of the present investigation and the 
summary results from Table 6 Indicate that the sources for the majority of 
the TID events are located in polar region. Some type of auroral activity 
occurs almost every night. Magnetic storms as well as particle precipitation 
cause intense auroral activity. The most Important mechanisms for generating 
atmospheric gravity waves operate In the auroral zone; neamely, joule heating 
due to auroral electrojets, heating due to particle precipitation and the 
motions Imparted to the neutrals by the motions of the charged particles. 

These source mechanisms may be responsible for generating large scale TID*s 
which can t>ropagate long distances. Even if there are no magnetic stoma 
to generate currents and other disturbances, still the night time auroral 
zone Is active with visual aurora forms moving with supersonic speeds, which 
appear to be responsible for auroral Infrasonlc waves (Wilson 1967). Particle 
precipitation Is a regular phenomena which in turn causes heating and move- 
ment of neutrals by the charged particle motions. Several cumulative effects 
operating In the auroral zone during disturbed and quiet conditions could be 
responsible for generating waves. Active auroral conditions during the night 
apparently cause higher speeds and longer wavelengths for night time TID 
events. 
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VI. CONCLUSIONS 


The following are some of the conclusions drawn from the study. 

1. During summer and winter seasons In both hemispheres, the TIDs 
travel away from the winter polar regions. During equinoctial months, 
disturbances in the northern hemisphere travel approximately westward while 
those In the southern hemisphere travel eastward. 

2. The daytime speeds are lower and periods are usually shorter 
than night time values. The daytime values are in the range of 70- 140 m 
sec~* and the night time values are In the range of 140- 340 m sec"^. The 
TID activity appears to be predominant during the winter season. 

3. The large scale TIDs are usually either accompanied by a 
magnetic storm or associated with some form of auroral zone activity. The 
medium scale TID’s are not generally correlated with either auroral or 
magnetic activity. 

4. The observed seasonal variation In TID propagation direction may 
be explained by mesospheric temperature variations. Observational data 
(Champion 1957) Indicate that there is an appreciable weakening of the 
winter duct around 80 km due to mesospheric warnings. These temperature 
variations In the mesosphere cause leaky duct conditions which allow the 
gravity wave energy to penetrate to F-reglon levels. 

5. Auroral zone activity appears to be the reason for the longer 
periods and higher speeds of the night time TIDs. 
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VII. SUGGESTIONS AND PROPOSED STIT)IES 

Evidence for relationships between meteorological and Ionospheric 
phenomena has been advanced by several Investigators (for example. Baker 
and Davies, 1968; Georges, 1968). A CW doppler sounder array was In 
operation at Huntsville, Alabama, during 1967-68 to investigate the coupling 
between the troposphere and the Ionosphere at times of static test firings of 
Saturn rocket engines to see whether the acoustic energy released at the 
ground level reached ionospheric levels. The static test firings of Saturn 
rocket engines did not provide any conclusive evidence of acoustic energy 
leaking from the ground to the Ionospheric levels (W. T. Roberts, Personal 
Communication). However, the CW doppler data (provided by W, T. Roberts) 
did Indicate an evidence of coupling between tropospheric source and the 
Ionosphere at the times of severe local thunderstorm activity. About 
80% of high frequency CW doppler fluctuations with periods less than or 
equal to six minutes, were well correlated with the local thunderstorm 
activity. Figure 13 Is a typical event showing the thunderstorm correlated 
doppler fluctuations. 

The present three dimensional array of the NASA/MSFC Is more 
sophisticated and better designed to study the coupling between tropospheric 
sources, le. thunderstorm effects and the Ionosphere. The North Alabama 
region la generally very active for severe thunderstorms and tornadoes and 
our array Is Ideally located to study the statistical nature of this apparent 
correlation. 


Pa 
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The atmospheric layers between 90-200 km (where the temperature begins to 
climb and the molecular weight changes)ls rather poorly understood. These 
changes produce a zone where acoustic cut-off frequency (W^) can become 
less than Brunt-Vaisala frequency (W|,). This is so called anomalous zone 
where > Wg (Johnston 1967). The dispersion diagrams (Tolstoy 1963) 
show that there exists a region where one cannot distinguish between 
acoustic and gravity waves. The dispersion characteristics of the naturally 
occurring wave trains (TlDs) whose propagation characteristics are controlled 
by the atmospheric structure between 90-200 km may be used to investigate 
thermospheric wave propagation. 

The phase height fluctuations in the ionosphere may be investigated 
using the CW doppler data and power spectral analysis. The variation of 
spectral peaks of the phase height fluctuations provides ameasure of the 
Burnt-Valsala frequency near the radio reflection height. The statistical 
properties of the CW doppler sounding at variable transmitting frequencies 
can be used to monitor certain physical parameters of the upper mesosphere 
and lower thermosphere. 

I 

There are no three dimensional observations of TIDs which provide 
information on wave front tilts and true phase velocities. The NA.SA/MSFC 
doppler sounder array data on the three dimensional TIDs and may be used 
to investigate wave characteristics seasonal and diurnal variations of the 
travelling ionospheric disturbances. 
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Figure 1 Graphical representation of wave front moving 
across the doppler sounder array. 






Figure 2 Graphical representation for computing phase velocity 
of the wave front. 
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Figure 4a. Typical TID observed on Sept. 23, 
1973, on 4.080 MHz operating 
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Figure 4b. Typical TID observed on Sept. 23 
1973, on 4.759 MHz operating 
frequency. 
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Figure 4c. Typical TID observed on Sept. 23, 
1973, on 5.734 MHz operating 
frequency. 
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Typical TID observed on Oct. 1 
1973, on 4.080 MHz operating 
frequency. 
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Figure 7a. Typical TID observed on Oct. 2, 
1973, on 4.080 MHz operating 
frequency . 
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Figure 8a. Typical TID observed on Oct. 3, 
1973, on 4.080 MHz operating 
frequency. 
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Figure 11. Diunal variation of TID 

direction of propagation in 
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Figure 13. CW Doppler record showing the 

thunder storm correlated doppler 
fluctuations. 



Figure 14. Block diagram of transmitting 
system. 
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Figure 15. Block diagram of Oscillator 
chasls. 
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Figure 16. Block diagram of receiving 
system. 
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TABLE 1 

?KAS£>PATR SOUNDER EQUIPMENT CHARACTERISTICS 


1* Transmitter Power Output: 2.5 kw peak emitted power 

2* Prequancy Range: 2 to 26 MHz, stepped manually or programmed 

(2 channels only) 

3. Programmed Frequency Steps: 

Selectable: 1, 10, 20 40 or 100 KHz steps per prf, 2, 4, 8 
or 16 prf. 

4. Pulse Repetition Rates: 20 or 40 per second per channel. Maximum 

of 4 channels available (Sequential pulsing 
Is used). 

5. Pulse Width: Selectable: 100, 200, 300 or 500 microseconds 

Selectable: Square or Gaussian pulse 

6. Recorded Data: Digitally processed for magnetic tape, two programmed 

receiver channels of: 

a. Pulse group delay (1 microsecond resolution) 14 bits (t) 

b. Pulse phase delay (3.6 degree resolution) 7 bits (4>) 

c. End of block frequency 16 bits each channel 

d. Julian day 16 bits 

e. Time of day (hour-minute- second) 24 bits 

f. Parity check for Items (1) and (2) above 

g. Data recorded in blocks of 600 pulse groups of group delay, 
phase delay, frequency, day and time. Each block consists 
of 4,813 bytes (4,812 6-bit bytes equals 3,609 8~bit words 
pulse one 6-bit longitudinal parity check word). 

7. Processed real-time delay: Two programmed receiver output channels 

of: 

a. pulse group delay (range to 900 km .... 6 milliseconds) 

b. pulse group differential with time 

c. pulse phase (+ 180 degrees relative phase) 

d. pulse phase differential with time 

Four receiver output channels of : 

a. pulse group delay, or 

b. pulse phase (+ 180 degree relative phase), or 

c. swept frequency lonograms (two channels) 



TABLE 2 


FORMULAS CONNECTING THE FREQUENCY VARIATIONS 


Equation and Comment 

Aaaumptlons 

e dt 

Baalc Equation 

Time and Space Variation of ionosphere 
alow compared to f (SVI) 

Af.l iMds 
cf \i at 

1 

SVI 

Concentric ionosphere (Cl) 
No Magnetic Field (NF) 

No Collisions (NC) 

Af • - 2fVv . ^ 
c 

Specular Reflection 

SVI, NF 

Stationary Transmitter and Receiver 
(STR) 

Vertlcal-Incldence Propagation (VI) 

A £ « 

Nondevlatlve layer 
(v « 1) 

SVI, NF, STR, NF, NC, Cl 

A f « AH 
Hall Drift Model 

Current flowing in the E-reglon 
produces fluctuating magnetic field at 
surface of Earth 

i( * 


SuddAn Commencement Model 

i 

A f « AHf" 

Alfven Wave Model 


. 1 h 

2v* dt* 
n 


. 1 i 

2v* f dt 
h 







TABLE 2 (continued) 

FORMULAS CONNECTING THE FREQUENCY VARIATIONS 


t * Time 

f • Transmitter frequency 

P * li cos a da phase path 

a " Angle between ray path and wave normal 

s ■ Ray~path lengthy measured from transmitter to point of Interest 
c " Speed of light, 3 x 10® m/s 

N “ Electron number density 

y “ Refractive Index given by the Apple ton-Har tree equation 

V " Velocity of reflecting surface. Subscripts v and h represent 
vertical and horizontal velocities, respectively 

AH ■ Fluctuation in the horizontal component of geomagnetic field 

h ■ Height of reflecting surface 




TABLE 

DISPEBSIOfi FOmiAS POE GEAVm HAVES AHD VHISTLEES 


Be£rmctlv« Index 

Uniting refrectlve index 
(X » Y 8ln it » 1) 

Critical prepay t Ion angle 

Angle a between wave nonal and 
ray (for long waves) 


Gravity Wave Foraulas 


Whlntler Have Fonnilaa 


- X - 1 


aln^ ♦ - I 


.* • 


w * 8la* t 


t* 8la* ^ 
“V In “ 


ao*) 


<IU) 


U28) 


j* - 


i* = 


|T cos 6|-1 


% 


{T cos 6( **1*^ 

C08 6 _ •in 


tan a • - cot 


u - C — sin 
“*a 


(13a) tan a tt 1/2 tan B 




“H 


(10b) 

(Ub) 

(I2b) 

(13b) 


Group velocity for very long waves 


(14a) 


o ■ 2c 


e 


(14b) 



OCHPAUSOB BBTBESr AVfOSFHEKlC BATES AHQ MACHETDIOHIC EADIO BATES (XSOIHEBMAL i^SHEBE) 


Aeo^tlc VoTca (IsotlwEMl ataosphcre) 




Radio Waves 


Aooafttlc cutcff frequency 

1 

3 


(U) 

Plasna frecpieiicy 

M * - 5sL 
ecP 

(lb) 

firunt-Valsalft freqiieiicy 

g H 


(2«) 

Oyrofrequency 


(2b) 


c. ^ 

X • 

ii) 


(3s) 


“u 1 

» - 

(3b) 


“g 

T - 


(4s> 


r . 

(*b) 

tefractiv* Index vitb 


X 

(5s) 

Refractive index with 

y* ^ cV.j _ 2X (l_- X) 




(5b) 

Srevlty 

M* 1 - T‘ 

sla* ^ 

aagnetlc field 

u* 2(1>Z)-T*sln*8t [T*sla*S+4T*(l-X)*««s*ei*^* 

Refractive Index vlth 

W » 

W* - 1 - I 


(«s) 

Refractove Index with 

U » Ug 

M 

1 

rH 

1 

#8 

(6b) 

Vertical propagatlcm 
(X < 1) 

U* - 1 - I 


(7s) 

Transverse propagation 
(ordinary wave) 

M* - 1 - X 

(7b)' 

Angle a betveen vave 


Jll 

(8s) 

Angle a between vave 

tjm a- - - sin 0eo«9 

(8b) 

vector k and the ray 

*“ * t* sin* *-l 


vector k and the ray 

[t* Bia"e+4T*a-X)* cos*ei>^* 

Group velocity (u) with 

y « 1 

- C* (1 - X) 


(9s) 

Croup velocity with 
t « 1 

“ c* (1 - X) 

(9b) 



TABLE 4 


Investigator 

Location 

Technique 

Year 

Observational Characteristics 

Suggested Source and Cotments 




NORTHERN HEHISPEERE 


Llszka and 
Taylor (1965) 

Irons 

69.6°N, 

19®E 

Faraday 
rotation of 
signals 

1963 

200-300 m. sec.~\ travel 
southward 

Auroral zone source. Corre- 
lation with magnetic activity. 

Goodwla (1968) 

Baker Lake 

64.3°N, 

96°W 

lonosonde 

Oct-Nov 

1962 

Horizontal and vertical wave- 
lengths range 800 - 13 km and 
travel away from winter noles. 

No correlation with geomagnetic 
activity. Associated with D- 
region absorption. 

Rlostermeyer 

(1969) 

Llndau 

51.7°N, 

10°E 

Spaced 

lonosonde 

Winter 

1960- 

1961 


Tropospheric Source 

Thomas (1959) 

Swansea 

51.6°N 

4°tf 

Spaced 

lonosonde 

June 

1951- 

Aug. 

1953 

40 - 280 m. sec."*travel SE 

TlD’s travelled with no change 
in characteristics over a 
distance of 650 km. Velocity 
increased with geomagnetic 
activity. 

Bramley and 
Ross (1951) 
Bramley (1953) 

Slough 
51.5°N, 
0,5° E 

Pulsed phase 
comparison 

1948- 

1952 

Daytime observations 35-339 m. 
sed”** 50-400 km horizontal 
wavelength 

The F-region echo showed quasi- 
period ic fluctuation of 3-30 
min. 

Vassuer and 

Waldteufel 

(1969) 

St. Santln- 

Nancy 

48.7°N, 

6°E 

Thompson 

Scatter 

Sept. 

1967 

Vertical wavelength 520 km, 
horizontal wavelength 
4000 km 

Auroral zone source. Correlates 
vrith magnetic activity. 


TABLE 4 (continued) 


Investigator 

Location 

Technique 

Tear 

Observational Characteristics 

Suggested Source and Conments 

Castle and 
Faynot (1964) 

Garchy 

47.3°N 

3.5°E 

Top and 
bottom 
iono sondes 

April 

1963 

100-200 m. sec.“^ 
Irregularities move from 
higher to lower latitudes 
under Sq currents 

Source located in high 
latitudes 


Elkins and 
Slack (1969) 

Boston 

42.50N 

123®E 

Transmissions 
from geo- sta- 
tionary sat- 
ellites 

1967 

50-120 m. sec. ^ travel 
toward equator, horizontal 
wavelengths 50-100 km 

Injection of energetic 
electrons into auroral 
zone ionosphere. 

Davies and 
Jones (1971) 

Boulder 

40°N, 

105®W 

CW doppler 

sounder 

array 

1967- 

1968 

145 m. sec.”^ velocity 
vector directed away from 
winter poles 

Sources located in polar 
regions. 

Bunsucker and 
Tveten (1967) 

Boulder 

40°N, 

105°W 

H. F. Back- 
scatter radar 

Winter 

1963- 

1964 

80-325 m. sec. ^ travel 
away from poles 

52% of the events preceeded 
by auroral activity. 

Davis and 
da Rosa (1969) 

Palo Alto 
37.30 N, 

122 ®W 

Faraday 

rotation 

from 

satellites 

Feb. 

1967 
Nov . 

1968 

440-600 m, sec,”^ travelling 
N-S direction horizontal wave- 
lengths 1500 km 

Source auroral oval during 
polar substorms. 

Chan and 
Vlllard (1962, 
1964) 

Palo Alto 

37 . 3 ON, 

I 22 OW 

CW doppler 
sounders 

1960-61 

400-765 m. sec.“* 30-50 min, 
periods horizontal wavelengths 
1300-1800 km 

TlD^s occurred following 
sudden connnenceinent of geo- 
magnetic storms. 

Georges 

(1968) 

Little Rock 

34 . 50 N, 

92°W 

CW doppler 

sounder 

array 

Aug, 

1966 
Feb. 

1967 

Medium scale < 300 m, sec,~\ 
large scale > 300 m, sec,’*' 

Large scale irregularities 
origin near poles 


TABLE 4 (continued) 


Invest leator 
Thome (1968) 

Location 

Technique 

Tear 

Observational Characteristics 

Suites ted Source and Comnents 

Puerto Rico 
18.3°N, 

60°W 

Thompson 

scatter 

Summer 

1964 

200-300 m. sec.”^ travel away 
from poles 

Auroral zone source 

Reddi and 
Rao (1971) 

Waltalr 

17.5°N, 

83°E 

Pulse phase- 
path sounder 

1966- 

1968 

40-230 m. sec. periods range 
3-40 min. 

Group speeds are twice the 
the phase speed. 





SOUTHERN HEMISPHERE 


King (1966) 

Hallett 

70.30s. 

170. 30 E 

lonosonde 

Summer 

1968 

200 m, sec.“^ travel away 
from pole 

Infrasonic waves and auroral 
zone source 

Titheridge 

(1968) 

Invercargill Faraday 
46.3®S, rotation 

168. 5 OE from satel- 

lites 

June- 

Nov, 

1962 

50 m, sec.”^ Winter 
periods 15-45 min . , 

Summer periods 15-70 min. 

Travelled away from south 
pole. Source located in 
polar region. 

Wright 

(1961) 

Christ- 
church 
43.5°S, 
172. 50 E 

lonosonde 

1957- 

1960 

406 - 690 m. sec." only NS 
component of the velocity 
was measured. 

TID's observed to a horizontal 
distance of 1100 tan. Positive 
correlation with magnetic 
activity , 

Munro 

(1958) 

Sydney 
34°S, 
151. 2°E 

Three fixed 
frequency 
pulsed trans- 
mitters 

April 

1948- 

March 

1957 

120 - 150 m. sec."' Travel 
away from winter poles 

No correlation with magnetic 
activity 

Heisler 

(1963) 

Sydney 
34°S, 
151, 2°E 

Spaced 

lonosondes 

1952- 

1954 

97 - 207 m. sec,"' travel away 
from winter poles, 

TID’s travelled with no change 
in characteristics over a 
distance of 3000 km. No 
variation of speed with height 


TABLE 4 (continued) 


InvestiRator 

Baker and 

Gledhlll 

(1965) 

Location 

Grabms 

Town 

33.3°S 

26.5°E 

T8c1ini<Lue 

Tear 

Observational Characteristics 

Suggested Source and Comments 

lonosonde 

August 

1960 

116 m. sec.“^ westward 
direction 

Source may be in sunrise 
conditions of P-Region. 

Price (1953, 
1955) 

Perth 

32^5, 

116°E 

Fixed fre- 
quency 
pulsed 
transmitters 

1949- 

1952 

90 - 330 m. sec.“' winter 
direction 0^ — 6® Elq , N . 

Summer direction 90° — 180° 

Eg. H. 

No correlation with 
auroral activity of K-lndex 

Bowman (1965, 
1968) 

Brisbane 
27.5°S, 
135. 5°E 

lonosonde 

Winter 

1958 

360 m. sec."^ horizontal wave- 
length 25 - 350 km. Periods 
5-60 min. 

Sources are conjugate high 
latitude locations, radio 
and optical auroras and 
magnetic activity. 


TABLE 5 


GEOGRAPHIC LOCATION OF THE CW DOPPLER SOUNDER ARRAY 


Station 

Latitude 

Longitude 

NASA-MSFC Alabama 

34° 39' N. 

00 

o^ 

o 

o 

Ft. McClellan, Alabama 

33° 44’ N. 

85° 48’ W 

TVA Muscle Shoals, Alabama 

34° 46' N. 

85° 38' W 

TVA Nclkajack Dam, Tennessee 

35° 01' N. 

87° 38' W 


o-Z? 



APPENDIX 


The NASA-MSFC CW doppler sounder array consists of nine CW high 
frequency field transmitters covering a surface area of about 150 km^ 
with radio receivers located In the NASA-MSFC, Huntsville, Alabama. The 
received signals from all nine transmitters are mixed with stable 
oscillators of the same nominal frequencies, but are offset several 
Hz lower. The beat frequencies are then recorded on a slowly moving 
(one Inch per minute) magnetic tape. For data analysis, the tapes are 
played back at higher speeds into an audio spectrum analyzer. The higher 
play back speed (for example 15 Inches per second) Increases all frequencies 
on the tape by a factor of 900, Into the range of about 1 to 5 kHz. The 
output of the spectrum analyzer Is connected to a facsimile recorder, which 
gives a strip chart record of the Doppler shifted frequency components vs. 
time. 

A block diagram of the equipment for one transmitting station Is shown 
in Figure 14. Three ovenlzed crystal oscillators are built Into one chassis , 
as shown In the block diagram of Figure 15. Since each of the three 
transmitters use the same nominal frequencies, the actual frequencies must 
be slightly different so that they can be distinguished from each other 
at the receiver. An adjustment on each oscillator module allows the 
frequency to be offset by about + 2 parts per 10®. The actual frequencies 
have been set as follows: 
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Nominal Frequency 

Actual Frequencies 

Unit 1 

Muscle Shoals 
Alabama 

Unit 2 

Nlkajack Dam 
Tennessee 

Unit 3 

Ft. McClellan 
Alabama 

4.0125 MHz 

4.0125 MHz 

4.0125 MHz 

4.0125 MHz 


+ 2Hz 

+ 3Hz 

+ 4Hz 

4.759 MHz 

4.759 MHz 

4.759 MHz 

4,759 MHz 


+ 2Hz 

+ 3Hz 

+ 4Hz 

5.734 MHz 

5.734 MHz 

5.734 MHz 

5.734 MHz 


+ 2Hz 

+ 3H* 

+ 4 Hz 


At the receiving station, the oscillators have been set exactly at the 
nominal frequencies, so that In the absence of Doppler variations, the 
offsets of the transmitter oscillators are equal to the beat frequencies. 

The stability of the oscillators is specified at + 1 part per lO’ per day 
and + 1 part per 10® over the ambient temperature range of 0 to 55*^ C. 

These variations will not be noticeable, but over a period of a year, the 
oacillators could drift by as much as 1 Hz. Probably, all the oscillators 
at both the transmitting and receiving stations will drift in the same 
direction and the effects will tend to cancel. 

The transmitters are modified versions of the Heathkit DK 60B. The 
transmitting antennas are halfwave dipoles. In order to maximize vertical 
radiation and minimize horizontal radiation, the height of the dipoles above 
ground should be between one eighth to one fourth wavelength. Since all 
three antennas at a site are In one array, a height of approximately 30 feet 
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has proven to be practical. To null the horizontal radiation in the 
direction of the receiver, the transmitting antennas at the sites are 
oriented so that the array is broadside to the receiving site and the 
ends of the antennas point to the receiver site (NASA-MSFC) . Ground 
waves are not expected to be a problem since the sites are farther than 
40 kn from the receivers. If a ground wave is received, it can be 
identified by its stable beat frequency on the doppler records. 

A block diagram of the receiving and recording system Is shown in 
Figure 16, The reflected ionospheric echoes from all the nine field 
transmitters are received by three Haraniarlund SP-600-JX-17 communication 
receivers at three different nominal frequencies. These receivers use 
single conversion for channel frequencies below 7.4 MHz, and double 
conversion above 7.4 MHz. In either case, the internal crystal HF 
oscillator is used Instead of the variable frequency oscillator. It is not 
necessary that the HF oscillator In the signal receiver has the stability 
of the reference receiver to be mixed with the signal from the reference 
oscillator. Thus any variations in frequency disappear when the signal 
Intermediate frequency (IF) is mixed with the reference IF. 

In order to obtain the accurate timing of ionospheric disturbances, 
a time code is recorded on channel 4 of the data tape. When the tape is 
played back for spectrum analysis, the channel 4 output is mixed with 
whichever data channel is being played back, so that data and timing appear 
simultaneously on the records. 

The output of the timing signal generator is a sine wave during the 



first 20 seconds of each minute, and nothing during the last 40. This 
signal produces a dashed line on spectrograph records, with the leading 
edge of the dash making the beginning of the minute. The minutes are set 
Into groups of five by a variation In frequency of the output signal. The 
frequency Is 1.0 Hz for five minutes, then 7/8 Hz for 5 minutes, then 1 Hz 
for five minutes, and so on. This way the minutes can easily be counted 
from some reference time. To mark the beginning of each hour, a group of 
five minute marks which would ordinarily be 7/8 Hz Is lowered to 3/4 Hz. 
Thus, the frequency Is 3/4 Hz for the first five minutes, 1.0 Hz for the 
second five, 7/8 Hz for the third five, ..., 1.0 for the last five. For 
identifying the hour, the frequency of the signal during the first five 
minutes of the 12th hour (GKI) of the day Is lowered to 1/2 Hz Instead 
of 3/4 Hz. The algorithm for selecting the frequency is as follows: 

1.0 Hz when minutes digit ■ 5, 6, 7, 8, or 9 

7/8 Hz when minutes digit ■■ 0, 1, 2, 3, or 4 and 

tens minutes digit + 0. 

3/4 Hz when minutes digit “ 0, 1, 2, 3, or 4 

and tens minutes digit ® 0 
and hours digit + 12. 

1/2 Hz when minutes digit ■*0, 1, 2, 3, or 4 

and ten minutes digit - 0 
and hours digit •• 12. 


The tape recorders used to record data are designed for recording 



at one Inch per minute only, and will not play back or rewind. The capstar 
twlna at 2 RPM and is machined to a diameter of 1/2 n (0.1596 In). Channels 
1, 2, and 3 are used for data, with a timing signal recorded on channel 4. 
The data tapes are played back on a four channel tape recorder at a speed 
of 15 inches per second. Each data tape is played back three times; with 
each playback, one data channel is mixed with the timing channel. The 
channels are recorded as follows: 

Channel 1: 4.0125 MHz, low frequency 

Channel 2: 4.759 MHz, middle frequency 

Channel 3: 5.734 MHz, high frequency 

Channel 4: Timing signal 

The output from the tape recorder is spectrum analyzed using a 
Federal Scientific UA-6®, Ublguitous audio spectrum analyzer. The output 
from the spectrum analyzer is fed to a Gifft GPR facsimile recorder, 
which gives a record of frequency vs. time. The receiving and recording 
systexos are shown in Figure 17 . 



